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ABSTRACT 

We examine the radial distributions of stellar populations in the globular cluster (GC) M15, using HST/WFC3 
photometry of red giants in the nitrogen-sensitive F343N-F555W color. Surprisingly, we find that giants with 
“primordial” composition (i.e., N abundances similar to those in field stars) are the most centrally concentrated 
within the WFC3 field. We then combine our WFC3 data with SDSS u,g photometry and find that the trend 
reverses for radii > 1' (3 pc) where the ratio of primordial to N-enhanced giants increases outwards, as already 
found by Lardo et al. The ratio of primordial to enriched stars thus has a U-shaped dependency on radius 
with a minimum near the half-light radius. N-body simulations show that mass segregation might produce 
a trend resembling the observed one, but only if the N-enhanced giants are ~ 0.25 M G less massive than the 
primordial giants, which requires extreme He enhancement (T > 0.40). However, such a large difference in 
Y is incompatible with the negligible optical color differences between primordial and enriched giants which 
suggest AT < 0.03 and thus a difference in turn-off mass of AM < 0.04 M Q between the different populations. 
The radial trends in M15 are thus unlikely to be of dynamical origin and presumably reflect initial conditions, 
a result that challenges all current GC formation scenarios. We note that population gradients in the central 
regions of GCs remain poorly investigated and may show a more diverse behavior than hitherto thought. 
Keywords: globular clusters: individual (M15) — stars: abundances — stars: Hertzsprung-Russell and C-M 
diagrams 


1. INTRODUCTION 

There is strong evidence, both from photometry and spec¬ 
troscopy, that globular clusters (GCs) have large internal star- 
to-star variations in the abundances of light elements. While 
some stars in GCs display the same chemical abundance pat¬ 
terns observed in metal-poor field stars, a large fraction of the 
GC stars (often the majority) exhibit combinations of light 
element abundances that are unique to GCs (apart from a 
small fraction of halo stars that may have escaped from GCs; 
|Martell et al.||20lT| ). This includes enhanced abundances of 
He, N, Na, and Al and depleted abundances of C, O, and Mg. 


email: s.larsen@astro.ru.nl 

1 Based on observations made with the NASA/ESA Hubble Space Tele¬ 
scope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations are associated with pro¬ 
gram #13295. 


Large spreads in the abundances of heavier elements (such 
as Ca and Fe) are relatively rare, although a significant frac¬ 
tion of the GC population may exhibit small (but detectable) 
spreads in iron abundance at the level of ~ 0.05 dex (Carretta 
|et al.|2009a||Willman & Strader|2012] >. 

The large and correlated spreads in light element abun¬ 
dances point to proton-capture nucleosynthesis at high tem¬ 
peratures as the main source of the observed abundanc e 
anomalies ( |Cottrell & Da Costa||198l] |Langer et al.|| 1993] >. 
However, the site where the processing takes place, as well as 
the mechanism by which processed material is subsequently 
incorporated into new generations of stars, remain uncer¬ 
tain. The main candidates for the production site are mas¬ 
sive asymptotic giant branch (AGB) stars, where the rele¬ 
vant nuclear reactions take place at the bottom of the con¬ 
vective envelope during hot bottom burning (Ventura et al. 


2001[ D’Antona & Ventura| j2007[ ), or massive (single or bi- 
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nary) main sequence stars (Wallerstein et al. |1987 Brown 


& Wallerstein 

1993; 

Prantzos & Charbonner|2006; Decressin 

et al. 2007b c 

e Min 

c et akT2009). In the AGB scenario, the 


polluted material is lost from the surface of the stars via slow 
winds that remain trapped within the gravitational potential 
of the cluster. In order to explain the observed continuous an¬ 
ticorrelation between [Na/Fe] and [O/Fe], some fraction of 
the polluted wind material must be diluted with “pristine” 
gas, i.e., gas with th e sam e co mposit ion as the original (first- 
generation) stars (D’Antona & Vcntur aj2007| . In the “winds 
of fast rotating main sequence stars” (WFRMS) scenario, it 
is assumed that massive stars within GCs rotate near break¬ 
up speed. Processed material is brought to the surface by 
rotational mixing, lost via a mechanical wind, and then ac¬ 
cumulates in a disk around the star, where the second gener- 
ation of (low-mass) stars i s assumed to form (Decressi n et al.| 
2007a| Krause et al.|2013|>. Alternatively, interacting or merg- 


rng massive binary stars may provide an efficient way to lose 
large amounts of processed material that could be incorpo¬ 
rated into an enriched population (|de Mink et al.|2009 !. This 


last scenario is attractive because a large fraction of massive 
stars are indeed observed to be members of binaries that will 


interact during their lifetime (Sana et al.|2013]l. 

In L ^ fL r-L^. A G 1 D r,v^/4 \UI7D\zlC r.r.1 


In both the AGB and WFRMS scenarios, only a small per¬ 
centage (~ 5%) of the initial mass of the first generation is 
returned in the form of polluted material, which leads to a 
“mass budget” problem. The observed large fractions of pol¬ 
luted stars in GCs are accommodated by assuming that most 
of the first-generation stars have been preferentially lost, im¬ 
plying that GCs we re initially a factor of > 10 mo r e massive 
than they ar e now (Decressin et al. 2007a 2010[ Vesperini 


et al. |2010|i. Such a coprous mass loss is, however, diffi- 


cult to accommodate for the metal-poor GCs in the Fornax, 
WLM, and IKN dwarf galaxies, among which at least the For¬ 
nax GCs show the same anomalies as Galactic GCs ( [Larsen) 
et al. |2014bj). In these dwarf galaxies, we have found that 
the GCs currently account for at least 20% of the metal-poor 
stars, which is difficult to reconcile with a loss of > 90% of 
the initial cluster mass ( [Larsen et al.|2012l|2014a) i. Other dif¬ 
ficulties with these scenarios are that no young massive star 
clusters with extended, on-going star formation have yet been 
found (despite extensive searches which have included clus- 
ters with masses similar to those expected for young GCs; 
Bastian et al.||2013a[ IBastian & Strader||2014[ |Cabrera-Ziri| 
et al. ||20i4 2015| >, and that the embedded phase las ts much 
shorter than expected from the WFRMS scenario (Bastian 

let al.|2014j >. _ _ 

Alternatively, |Bastian et al.|(|2013b|> have suggested that the 
polluted material lost from interacting binaries may be swept 
up by accretion disks around low-mass stars. This “early 
disk accretion” scenario might provide a solution to the mass- 
budget problem, because (1) large amounts of polluted ma¬ 
terial are available from the interacting binaries and (2) the 
material is accreted onto pre-existing low-mass stars. In this 
scenario, there is thus only a single “generation” of stars. The 
time scale for the accretion to take place may, however, re¬ 
quire some fine tuning to ensure that a sufficient amount of 
ejecta are accreted and mixed while the low-mass stars are still 
in the convective phase (D’Antona et al. 2014| >. In the remain¬ 
der of this paper, we will generally refer to the stars with field¬ 
like composition as “primordial” and those that have modi¬ 
fied light-element abundances as “enriched” when discussing 
our observations, and thereby avoid implying a particular se¬ 
quence of events. 


All of the above scenarios predict that the enriched stars 
should be located preferentially in the central regions of the 
clusters. In the AGB scenario, the first (primordial) gen¬ 
eration is envisioned to expand following expulsion of gas 
left-over from the initial burst of star formation, after which 
wind material accumulates in the center via a cooling flow 
and forms a more centrally concentrated enriched population 
( |D’Ercole et al.12008] ). In the WFRMS scenario, the enriched 
stars are also expected to form preferentially in the central re¬ 
gions, because they form in the vicinity of mass-segregated 


massive stars (Decressin et al. 2008). A similar prediction 
is made by the early disk accretion model because accretion 
is more efficient in the central regions where the density is 
higher. These expectations appear to be borne out by obser¬ 
vations of some GCs, where a number of studies have found 


ter ( 

Norris & Freeman 1979; 

Carretta et al.;2009b; Kravtsov 

et a 

1. 2010; Lardo et al. [201 

'jMilone et al. 2012b). How- 

ever 

, the differences in spatia 

distribution are expected to be 


eventually erased by dynamical evolution. This should hap¬ 
pen first in the central regions of clusters, where the relaxation 
time is shortest ( Vesperini et al. 2013] >. Indeed, it has recently 
been found that the two populations in the cluster NGC 6362 
do not exhibit any differences in their radial distributions (jDa-| 
|lessandro et al.|2014| . 

To a large extent, then, the spatial distributions of stellar 
populations within most GCs observed to date appear con¬ 
sistent with theoretical expectations that are common to all 
formation scenarios. However, there may be more subtle dif¬ 
ferences between the predictions of different scenarios. One 
such difference concerns stars with intermediate composition, 
which are supposed to have formed out of diluted wind mate¬ 
rial in the AGB model and would therefore have been the last 
to have formed. In the early disk accretion scenario, the inter¬ 
mediate population would instead correspond to stars that did 
not pass through the densest part of the cluster. It is, there¬ 
fore, a clear prediction of this scenario that such stars should 
have an intermediate degree of central concentration. In the 
other scenarios it is less clear what to expect, but it seems 
plausible that the intermediate-composition population may 
be expected to be the most centrally concentrated in the AGB 
scenario, since it is the last to form. While this remains some¬ 
what speculative, it does suggest that interesting constraints 
on formation scenarios may be obtained by studying the ra¬ 
dial distributions in more detail. 

In light of the scenarios outlined above, we have used Hub¬ 
ble Space Telescope (HST) observations to examine the spa¬ 
tial distributions of red giants in the globular cluster M15 
(NGC 7078) as a function of their chemical composition. 
With a metallicity of [Fe/H] = -2.3 (Carretta et al. 2009a| ), 
M15 is one of the most metal-poor GCs in the Milky Way. 
While the internal spread inJEe/H] is small, estimated at 
o"rFe/Hi ~ 0.05 dex ([Carretta et al.|2009a[|Willman & Strader 
2012), M15 is similar to other Galactic GCs in showing large 


internal abundance variations of the light elements. Observa¬ 
tions of red giants have revealed the well-known Na/O anti¬ 
correlation, as well as a clear Mg/Al anti-correlation (Sneden 
et al. 1997). The abundances of C and N also exhibit large and 


anti-correlated variations with [C/Fe] an d [N/Fe] ratios vary- 
ing by ~ 1 dex and ~ 2 dex, re spect ively (Trefzger et al. 1983 


Cohe n et al.||2005] |Pancino et al.||20T0 ). In addition to the 


light-element abundance variations, M15 is one of a few GCs 
that are known to exhibit a spread in the abundances of the 
heavy (u-capture) elements Sr, Y, Zr, Ba, La, and Eu, which 
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however do not appear to correlate with the abundance vari¬ 
ations of the light elements and may be of a different origin 
( Sneden et al.|1997j|Otsuki et al.|2006[ |Sobeck et al.|2011 1 . 
In spite of relatively short exposure times, our data are sensi¬ 
tive to N abundance variations for stars at the base of the red 
giant branch (RGB). Combined with the high luminosity of 
the cluster (My ~ —9.1; Harri s] 1996[ >. this yields a sample of 
more than 1300 red giants, which allows us to examine radial 
trends of the sub-populations in some detail. In this paper, we 
report the (unexpected) results of our investigation. 

2. OBSERVATIONS AND DATA REDUCTION 

As part of our HST program to study the globular clusters in 
the Fornax dwarf spheroidal galaxy (|Larsen et al.[2014b here¬ 


after Paper I), we obtained short exposures of M15 in the same 
filters as those used for our main program: F343N, F555W, 
and F814W (Program ID: GO-13295, P.I.: S. S. Larsen). 
These observations exploit the well-established sensitivity of 


(Hesser et al.j 1977 

Grundahl et al] 

2002; 

Yong et akl|2008 

Sbordone et al!|20 

Monelli et al. 

2013 

). The integration 

555Wand F814W 

time was 2x350 s in 

F343N and 2xlC 

s in F 


with M15 centered on CCD #2 of the Wide Field Camera 3 
(WFC3). The post-flash option was used to mitigate the effect 
of charge transfer losses by increasing the background level 
to 10 counts per pixel. All exposures were obtained within a 
single orbit. Note that these exposures were not designed to 
be very deep, but were only intended to reach stars on the red 
giant branch. 

The pipeline-reduced images were corrected for 
charge transfer inefficiencies with the program 
wfc3uv_cterevers£] We then used the astrodrizzle 
code to align, combine, and resample the two exposures 
in each filter to a uniform pixel scale of 07040 per pixel. 
Point-spread function (PSF) fitting photometry was carried 
out with ALLFRAME ( |Stetson| 1994[ > and calibrated to standard 
STMAG magnitudes as described in Paper I. 

In addition to our own WFC3 data, we use imaging of 
M15 in the F606W and F814W filters obtained with the Ad¬ 
vanced Camera for Surveys (ACS) on HST as part of the ACS 
Galactic Globular Cluster Survey (ACSGCS; Sarajedini et al. 
|2007[ >. The ACSGCS data consist of short exposures with 
integration times comparable to those of our F555W/F814W 
data (15 s in each filter), as well as deeper exposures (4 x 130 
s in F606W, 4 x 150 s in F814W) that allow accurate photom¬ 
etry for stars well below the main sequence turn-off. We did 
not carry out photometry on these images ourselves, but use 
the catalogs published by the ACSGCS team ((Anderson et al.) 
[2008). 

Throughout thi s pap er we assume a distance of 10.3 kpc 
(van den Bosch et akj|2006jl, along with a foreground ex¬ 
tinction of Ay — 0.30 m ag (jSchlafly & Finkbeiner 2011b . 
Using the Cardelli et al. ( 1989j )' extinction law, this yields 
Afm 3 n = 0.483 mag, Afsssvy = 0.312 mag, Ap^w = 0.276 
mag, and Apsiw = 0.178 mag in the HST filters. 

3. RESULTS 

3.1. Photometric evidence for a spread in the N abundance 

Figure [T] shows the (F343N-F555W, F555W) color- 
magnitude diagram for the lower part of the M15 RGB (a 
larger color/magnitude range, as well as the (F555W-F814W, 

2 http://www.stsci.edu/hst/wfc3/tools/cte_tools 



Figure 1. (F343N-F555W,F555W) color-magnitude diagram showing the 
lower RGB of Ml5. Symbols are color-coded according to the division into 
normal, intermediate, and strongly enhanced N abundances (group A, B, and 
C). The thick black line is a median ridge line. 

F555W) CMD, are shown in Paper I). Here, and in the follow¬ 
ing, we only include stars in the radial range 4 " < R < 130" 
for which the ALLFRAME photometry has^-J < 3 in the F555W 
images. Closer to the center, crowding prevents accurate pho¬ 
tometry and for R > 130" there are very few stars although the 
outermost corner of the WFC3 mosaic is nominally located at 
R ~ 150" from the center of M15. Table [T1 lists the photom¬ 
etry for all stars brighter than F555W = T9 (Mp^ssw = 3.6) 
that meet these criteria. 

As in Paper I, we will generally exclude stars brighter than 
■^F 555 w = 1 from our analysis, because their surface abun¬ 
dances may have been modified by deep mixing ((Gratton et al. 
|2000| l. However, because of the better S/N of the M15 data, 
we can obtain good photometry for somewhat fainter RGB 
stars (1.0 < Mp 555 w < 3.0) than in Paper I, where we adopted 
a limit of Mpsssw = 2.5 mag. The symbols for stars in this 
magnitude range are color-coded in Fig. |T| according to their 
offset in the (F343N-F555W) direction, as discussed further 


below (Sec. 3.4 1 . On first inspection, we note that the spread 


of the RGB stars in F343N-F555W is far greater th an th e pho¬ 
tometric uncertainties, which are ~ 0.02 mag (Sec. |3.3[ ). Fur¬ 
thermore, the spread does not depend significantly on the lo¬ 
cation along the RGB, again consistent with most of it being 
real. 

We also plot model colors for N-normal c ompositi on 
and the N-enhanced 

[20TT ) 


‘CNONal” mixture (Sbordone et al7| 
A([C/Fe], [N/Fe], [O/Fe], 


which has A([C/Fe], [N/Fe], [O/Fe], [Na/Fe]) = 
(-0.6, +1.8, -0.8, +0.8) dex relative to standard (cr-enhanced) 
composition. The colors were computed for a 13 Gyr 
isochrone ( [Potter et al.||2007] > with [Fe/H] = -2.3 and 
[q/ Fe] — +0.4 by integrating ATLAS 12/ SYNTHE model spec¬ 
tra ( |Sbordone et al.|2004[ Kuruc z|20051 > over the filter trans¬ 
mission curves. The two model lines are clearly separated, 
with the N-enhanced models being redder by about 0.16 mag, 
and the observed F343N-F555W colors span a range compa¬ 
rable to, or somewhat greater than, the separation between 
the models. This is in agreement with the spectroscopically 
measured N abundance spread in M15 of about 2 dex (jCo-| 
|hen et ak||2005j). We note that the models do not match the 
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Table 1 

ALLFRAME photometry of M15 


ID 

X 

Y 

R.A. 

Decl. 

R 

F343N 
mag err 

F555W 
mag err 

F814W 
mag err 

102 

658.22 

51.70 

322.508246 

12.154236 

70.86 

18.607 

0.036 

17.975 

0.017 

18.345 

0.018 

108 

534.34 

54.74 

322.509641 

12.154211 

74.80 

18.874 

0.037 

18.318 

0.020 

18.688 

0.019 

272 

304.00 

84.04 

322.512246 

12.154424 

81.94 

19.176 

0.043 

18.904 

0.031 

19.499 

0.026 

338 

1138.58 

91.72 

322.502863 

12.154904 

55.56 

18.820 

0.029 

18.430 

0.012 

18.887 

0.010 

353 

952.61 

93.77 

322.504955 

12.154838 

60.74 

18.518 

0.021 

17.922 

0.020 

18.274 

0.020 

355 

1363.93 

93.60 

322.500329 

12.155031 

49.90 

18.995 

0.016 

18.647 

0.020 

19.138 

0.019 

374 

1236.21 

95.60 

322.501766 

12.154993 

52.90 

19.285 

0.027 

18.981 

0.014 

19.604 

0.020 

389 

1194.67 

97.66 

322.502235 

12.154996 

53.90 

19.337 

0.025 

18.964 

0.011 

19.632 

0.018 

406 

1160.12 

100.44 

322.502624 

12.155010 

54.72 

19.111 

0.027 

18.748 

0.013 

19.263 

0.018 

463 

1605.34 

105.80 

322.497620 

12.155280 

44.70 

17.613 

0.023 

16.942 

0.016 

17.296 

0.012 


Note. — X and Y are the coordinates in the drizzled CCD frames. R is the projected distance from the cluster center 
(in arcsec). For each filter, the magnitude and error estimate by ALLFRAME are listed. (This table is available in its entirety 
in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content). 


sub-giant branch perfectly but appear slightly too bright/blue. 
A better match to the subgiants can be achieved by increasing 
the extinction correction by A A v ~ 0.02 mag and the assumed 
distance by -0.25 kpc, which is well within the 0.4 kpc uncer¬ 
tainty (van den Bosch et al. |2006| . However, no conclusions 
in this paper are affected by these small adjustments and in 
what follows we simply keep the literature values. 

Apart from the color variations in F343N-F555W that 
arise from variations in the light element mixture, the color- 
magnitude diagram can also be affected by He abundance 
variations ( |Norris|2004||Salaris et al.||2006) l. While the vari¬ 
ations in F343N-F555W color are mainly an atmospheric ef¬ 
fect, due to the molecular absorption bands in the UV (pre¬ 


dominantly the NH band near 3370 A), He abundance varia¬ 
tions modify the internal structure and effective temperature 
of the stars and therefore also affect optical colors. These 
effects are illustrated in Fig. [2] for the WFC3 filters used in 
our program. We have combined isochrones for Y = 0.25 
and Y = 0.40 with ATLAS12/SYNTHE model atmospheres 
and synthetic spectra computed specifically for these Y val¬ 
ues and normal and CNONal light-element mixture. A fixed 
iron abundance, relative to hydrogen, of [Fe/H] = -2.3 
was assumed in all cases. As found by other authors, the 
F555W - F814W colors of stars on the lower RGB are virtu¬ 
ally insensitive to the CNONa abundances, but become bluer 
for increased He abundance (upper panel). In contrast, the 
F343N - F555W colors are very sensitive to the CNONa 
abundances and become redder as the N abundance increases 
(lower panel). This shift is much greater than the shift to¬ 
ward blue colors caused by a high He fraction. We thus 
expect that a N-enriched population will indeed have red¬ 
der F343N - F555W colors than a population with primor¬ 
dial composition, even if the N-enriched population is also 
strongly He enhanced. 


3.2. Differential reddening 

In addition to the mean foreground reddening toward M15, 
reddening variations on smaller scales may be present within 
the held of view of the HST cameras and could potentially 
affect analyses of spatial trends. Milone et al.]( 2012c] > used 
the mean colors of main sequence stars to map the redden- 
ing across GCs and found significant variations across the 
HST/ACS held of view in regions of high foreground redden¬ 
ing (E(B - V) > 0.1). For GCs in regions of lower foreground 
reddening, small reddening variations are difficult to disentan¬ 
gle from systematic variations in the photometric zero-points 




Figure 2. Model colors for different He and C,N,0,Na abundances. The 
isochrones have an age of 13 Gyr and [Fe/H] = -2.3. The F555W-F814W 
colors are essentially independent of the CNONa mixture, but sensitive to 
He abundance, whereas the F343N-F555W colors are mainly sensitive to 
CNONa variations. 
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F606W-F814W 

FigureS. (F606W-F814W,F606W) CMD for the upper main sequence of 
Ml 5. The ridge line and the limits used for selection of stars for the reddening 
map are shown. The arrow is the reddening vector for Ay = 0.3 mag. 
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Figure 4. i?(F606W - F814W) reddening map based on F606W - F814W 
colors of main sequence stars. The r.m.s. variation in £(F606W - F814W) 
across the field is it = 0.011 mag, with peak-to-peak variations of ~ 0.05 
mag. 


across the field that may be caused by uncertainties in the PSF 
modeling. Indeed, Anderson et al. ( |2Q08] > found systematic 
variations of about +0.01 mag in the mean F606W-F814W 
colors of main sequence stars even for GCs in regions of very 
low foreground reddening. Since the reddening toward M15 
is not entirely negligible, it is worth examining whether there 
is evidence of differential reddening in the HST data. 

Our WFC3 data are not deep enough to allow reliable mea¬ 
surements of main sequence stars, so we used the ACSGCS 
photometry to map variations in the colors of main sequence 
stars. A small fraction of the WFC3 field falls outside the 
area covered by the ACSGCS data, but this only affects re¬ 
gions with a projected distance of > 140" from the center of 
M15. We first defined a ridge line by computing the median 
F606W - F814W colors of main sequence stars as a function 
of F606W magnitude in bins of 0.1 mag. We only included 
stars with a photometric error in F606W of less than 0.0 2 mag 
and with a “quality-of-fit” parameter (Anderson et al. 2008| 
of qfitV < 0.3. Figure [3] shows the upper part of the main 
sequence with the ridge line overplotted. 

We then calculated the offset A 606-814 in the F606W - 
F814W color with respect to the ridge line for each star in 
the range 4.0 < M F 606 w < 6.0. Stars with |A 6 O 6 - 814 I > 0.06 
mag were excluded (dashed curves in Fig. [3J, leaving a total 
of 29515 stars. The dispersion of these stars around the ridge 
line was cr = 0.020 mag in the F606W - F814W color. The 
A 606-814 offsets were converted to estimates of the reddening 
by taking into account the slope of the ridge line and the red¬ 
dening vector. 


£(F606W-F814W) = 


^606-814 



/ <5606-814 
\ <^606 


I Af(,Q6W 
\Af606W - A/T814W 



( 1 ) 


where is the (inverse) slope of the ridge line, evaluated 

locally at the magnitude of each star. This relation is valid 
for offsets that are small enough so that the curvature of the 
ridge line can be neglected locally. For any individual star, the 
offset from the ridge line is typically dominated by random 


photometric errors, but by averaging the offsets of many stars 
we could map systematic variations in the mean color across 
the field. 

The sky coordinates from the Anderson et al. |2008) cata¬ 
log were transformed to pixel coordinates in the WFC3 frame 
with the rd2xy task in the drizzlepac package. At each 
pixel in the WFC3 image, we computed a weighted aver¬ 
age of the reddening values of the surrounding stars, with 
weights given by a Gaussian function of the distance, w, = 

ex p(~5 ( 100 pixels)- ) for a star located n P ixe l s from a 8 iven 

point ( Larsen11996[ ). On average, this yielded about 50 stars 
per “resolution element” of the resulting map, although the 
stellar density obviously varies greatly across the field. 

The map is shown in Fig. [4] The color scale indicates the 
average computed reddening (£’(F606W - F814W)) at each 
position and the black points show the locations of the stars 
used to produce the map. The gap in the distribution of stars 
near the center of the lower CCD detector corresponds to the 
center of M15, and it can be seen that the upper right-hand 
corner of the WFC3 field falls outside the coverage of the AC¬ 
SGCS data. Note that the computed color may vary even in 
regions with few (or no) stars, depending on how the relative 
weights of the closest stars to any given point change with po¬ 
sition. This is seen in the upper right-hand corner, as well as 
in the region near the center. It is clear that the significance of 
any features in these empty regions is low. 

By comparing maps from two independent sub-samples, we 
found that the large-scale features in Fig. [4] were consistently 
reproduced: there is an overall gradient across the field with 
redder colors (that may be interpreted as higher reddening) 
in the bottom left part of the field and bluer colors in the 
top right-hand part. There is some evidence for a filamentary 
structure extending across the lower half of the field from left 
to right, whereas most structure on smaller scales appears not 
to be significant. The r.m.s. deviation in LfFfiOfiW - F814W) 
across the map is 0.011 mag, which is comparable to the vari¬ 
ations expected from uncertainties in the PSF modeling (An-| 
|derson et al.|2008 |Milone et al.|2012c] ). 

Apart from differential reddening and position-dependent 
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uncertainties in the photometric calibration, other factors that 
might affect the optical colors include He abundance, over¬ 
all metallicity, and binarity. If variations in these quantities 
were significant, we would expect them to show up as a ra¬ 
dial gradient with respect to the center of M15. Given that 
the structure in Fig. |4l is not obviously symmetric around 
Ml5, it seems unlikely that these factors contribute signifi¬ 
cantly to the structure in the map. However, it remains diffi¬ 
cult to disentangle reddening variations from variations in the 
photometric zero-points. The dispersion of the MS stars in 
^ 606-814 decreased from cr = 0.020 mag to cr = 0.017 mag 
when correcting the photometry for reddening according to 
the map, which is consistent with the 0.011 mag dispersion in 
E(F606W - F814W) across the field. However, this decrease 
in the dispersion of the corrected colors is expected whether 
the variations are caused by differential reddening or photo¬ 
metric zero-point variations. A decrease was also seen for the 
dispersion across the RGB; without any correction the disper¬ 
sion was cr = 0.020 mag in F606W-F814W, and when apply¬ 
ing the reddening map this decreased (slightly) to cr — 0.019 
mag. However, in the WFC3 data, the dispersion of the RGB 
stars in F555W-F814W actually increased (from cr = 0.024 
mag to cr = 0.026 mag) if a reddening correction was applied. 
This suggests that the color variations in Fig. [4] are primar¬ 
ily caused by instrumental effects (and thus do not reproduce 
between the ACS and WFC3 observations), rather than by dif¬ 
ferential reddening. 

Thus having considered the possible effect of differential 
reddening carefully, we have at the end chosen not to correct 
for it in our general analysis. However, when relevant we 
will comment on any differences that arise from including or 
omitting this correction. 

3.3. Artificial star tests 

To quantify the photometric errors in the WFC3 observa¬ 
tions, we carried out artificial star experiments following the 
general procedure described in Paper I. We started by generat¬ 
ing coordinate lists for a number of concentric annuli around 
M15. The annuli covered the radial intervals 100-200 pix¬ 
els, 200—400 pixels, 400-600 pixels, and 800-1000 pixels. In 
each annulus, 500 pseudo-random star coordinates were gen¬ 
erated by arranging the stars in a polar grid with a spacing of 
20 pixels in the radial direction and a spacing in the azimuthal 
direction that provided the desired total number of stars. A 
further random dither offset in the range -0.5 ... + 0.5 pixels 
was added to each coordinate in the x- and v-directions. The 
area of the innermost bin was too small to accommodate 500 
stars with a minimum separation of 20 pixels, so for this bin 
only 100 coordinates were defined. We then generated lists of 
F343N, F555W, and F814W magnitudes for the artificial stars 
by selecting the F555W magnitudes of the actual RGB stars 
in M15 and interpolating in the N-normal model isochrone to 
find the other magnitudes. The artificial stars were added to 
the images using the mksynth task in baolab ( jLarsen|1999| >, 
including a set of artificial PSF stars. The ALLFRAME pho¬ 
tometry procedure was repeated, and the artificial stars were 
recovered by requiring a match within a distance of 1 pixel 
from the input coordinates. This procedure was repeated four 
times. 

Figure [5] shows the CMD of the recovered artificial stars. 
By comparison with Fig. IT] it is clear that the dispersion in 
the artificial CMD is much smaller than the observed spread 
in the F343N-F555W colors. As in Paper I, we defined 
the A(F343N-F555W) parameter as the offset between the 



Figure 5. (F343N-F555W.F555W) CMD for synthetic stars populating a sin¬ 
gle isochrone. The broadening of the RGB due to photometric errors is far 
smaller than the observed width of the RGB seen in Fig.^j 

isochrone of N-normal composition and the observed F343N- 
F555W color at a given F555W magnitude. In Fig.|6]we com¬ 
pare the distributions of A(F343N-F555W) for the observa¬ 
tions of red giants in M15 with the artificial star tests for the 
magnitude range 1 < Mp 555 W < 3. Even though the radial 
density distribution of the artificial stars is not fully realistic, 
the histogram for the artificial stars is much narrower than for 
the observed RGB stars. Formally, the dispersion of the ob¬ 
served A(F343N-F555W) values is 0.059 mag, whereas the 
corresponding dispersion for the artificial stars is 0.018 mag. 
In A(F555W - F814W), the synthetic CMD has cr = 0.020 
mag, which is slightly less than the observed spread for the 
RGB stars (cr ~ 0.024 mag). This is consistent with an addi¬ 
tional ~ 0.01 mag variation from uncertainties in PSF model¬ 
ing and/or differential reddening. 

To quantify the dependency of the photometric errors on 
magnitude and radial position further, we carried out a sec¬ 
ond set of artificial star tests in which stars with fixed magni¬ 
tudes of F555W = 17, 18, and 18.5 (Mf 555 w ~ 1.6, 2.6, and 
3.1) were added to the images. The radial bins and numbers 
of stars in each bin were the same as previously described. 
The resulting dispersions in F343N - F555W (cnj_ v ) and in 
F555W - F814W (cr V -i) are listed in Table [2] for each radial 
bin and each input magnitude. As expected, crowding causes 
the errors to increase toward the center, especially for radii 
< 200 pixels (8"), but the errors remain much smaller than 
the observed color spread at all radii and magnitudes. We have 
included the recovered F343N-F555W distribution for the in¬ 
nermost annulus (4" - 8") and F555W = 18.5 as a thin (red) 
histogram in Fig. [6] Even for this “worst case” (where the ar¬ 
tificial stars are 0.1 mag fainter than our magnitude limit), the 
color distribution of the artificial stars is much narrower than 
the observed one. At all magnitudes and for all radial bins, 
more than 95% of the synthetic stars were recovered. 

3.4. Radial trends in the HST/WFC3 data 

In Fig. [7] we plot the observed A(F343N-F555W) versus 
projected distance from the cluster center. The visual im¬ 
pression from this figure is that the stars with the bluest 
A(F343N-F555W) colors (A(F343N - F555W) < 0.05) ap- 
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Table 2 

Photometric errors and completeness from artificial star tests 


Radius 

F555W = 

17 

F555W = 

18 

F555W = 18.5 


cry—I 

cru-v 

free 

CTy_i 

CTy-V 

free 

CTy_i 

cru-v 

free 

4"~8" 

0.013 

0.013 

0.99 

0.029 

0.032 

0.99 

0.045 

0.038 

0.99 

8"—16" 

0.011 

0.013 

0.99 

0.021 

0.024 

0.97 

0.025 

0.031 

0.98 

16"-24" 

0.011 

0.011 

0.99 

0.017 

0.018 

0.99 

0.021 

0.022 

0.99 

32"-40" 

0.009 

0.009 

1.00 

0.016 

0.015 

0.99 

0.019 

0.017 

1.00 


Note. — trv-i and cru_v denote the standard deviation of the recovered 
artificial star colors in F555W-F814W and F343N-F555W, respectively. / rec 
is the fraction of the input artificial stars recovered by the photometry proce¬ 
dure. 



Figure 6. Distributions of A(F343N-F555W) values for red giants in M15 
(filled histogram) and the synthetic CMD (unfilled black histogram). The 
thin red histogram shows the F343N-F555W distribution for artificial stars 
with F555W = 18.5 mag in the innermost annulus. A' 4” - 8" (cf. Table[2j. 
The red histogram has been shifted by 0.03 mag. 

pear to be the most centrally concentrated. There may 
also be a difference between stars with intermediate colors 
(0.05 < A(F343N - F555W) < 0.13) and the reddest (most 
N-enhanced) stars. 

Based on this first assessment of the data, we divided the 
stars into “N-normal”, “intermediate” and “extreme” sam¬ 
ples, corresponding to -0.03 < A(F343N - F555W) < 
0.05, 0.05 < A(F343N - F555W) < 0.13, and 0.13 < 
A(F343N - F555W) < 0.23, respectively. To avoid con¬ 
tamination by blue stragglers, we further required F555W - 
F814W > -0.6 (see Figure 4 in Paper I). These groups con¬ 
tain 32%, 50%, and 18% of the 1339 RGB stars in the mag¬ 
nitude range considered here and in the following we refer to 
them as groups A, B, and C. The different colors in Fig. [T] 
correspond to these three groups. The dashed, vertical lines 
in Fig. [7] show the median radial coordinates of each group; 
these are R me d,A = 23"3 ± 171, R me d,B = 2974 ± 171, and 
f^med.c = 3274 + 270 for the three groups, respectively (er¬ 
rors were estimated via bootstrapping). The group A stars are 
clearly more concentrated than the other groups. The artificial 
star tests showed that detection incompleteness is negligible 
over the magnitude- and radial ranges considered here. How¬ 
ever, the spatial completeness drops below 100% at distances 
of > 40" from the cluster center and is only 10% at our outer 


Radius [pc] 



Figure 7. The A(F343N-F555W) color offset (relative to the “Std composi¬ 
tion” isochrone in Fig. 17} as a function of projected distance from the center 
of M15. The vertical (red) dashed lines indicate the median radial coordinate 
of each color group. 

limit, 130". While this should not affect the comparison of 
the different groups in a strictly relative sense, the absolute 
values of the median radii are therefore not related to physical 
cluster properties in a simple way. 

Our A, B, and C groups are somewhat reminiscent of the 
primordial (P), intermediate (I) and extreme (E) populations 
defined by Carretta et a I.] ( j2009b[ > based on Na and O abun¬ 
dances. However, while Carretta et al. found a good cor¬ 
respondence between their spectroscopically defined popula¬ 
tions and CN-sensitive photometry, it is not clear that their 
groups are exactly equivalent to ours. Indeed, they find no “E” 
stars in M15. To avoid confusion, we therefore use a different 
naming scheme. We emphasize that the adopted division is 
not meant to imply the existence of three distinct populations. 
Indeed, there is no evidence for this in our CMD. The pho¬ 
tometry of Piotto et al. ( |2015} does suggest a bimodal struc¬ 
ture of the RGB in M15, but al so has a s ubstantial number of 
stars with intermediate colors. |Pancino et al.| ( |20T0) > found a 
bi-modal distribution of CN and CH band strengths from low- 
dispersion spectroscopy of main sequence stars in M15, while 
bimodality is less clear or absent in other data (jKayser et al. 
[20081 |Cohen et al.|2005HSneden et al.|1997] l. 

Figure [8j shows the cumulative radial distributions of the 
three groups. These confirm that the group A stars are the 
most centrally concentrated, followed by group B, and the 
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Figure 8. Cumulative radial distributions of the populations. The primordial 
(group A) stars are clearly the most centrally concentrated. 

group C stars are the least concentrated. From a K-S test, we 
find that the P values when comparing the radial distributions 
of the group A and B stars, the A and C stars, and the B and C 
stars are P A b = 1.4x10~ 3 ,Pac = 1.4x 10 -5 , and P B c =0.19, 
respectively. The A stars thus differ very significantly from 
both the group B and C stars, whereas the difference between 
the B and C stars is only marginally significant. This is con¬ 
sistent with the error estimates on the median radii above. 

These findings are robust to changes in the exact selection 
criteria and details of the analysis. In Table [3] we list the 
median radii and P-values for different magnitude cuts and 
other modifications to our analysis procedure. If the correc¬ 
tion for differential reddening is included, the differences be¬ 
come even more significant. We also tried excluding blended 
stars, here defined as stars that have a neighbor within the fit¬ 
ting radius (3 pixels) that is brighter than F555W + 2, where 
the F555W magnitude refers to the magnitude of the star it¬ 
self. This criterion removes about 5% of the stars and the 
P values increase somewhat. Because the blended stars are 
found preferentially in the inner regions of the cluster, the me¬ 
dian radii all increase slightly, but the overall differences are 
preserved. For brighter magnitude limits, the number of stars 
decreases and with it the statistical significance of the differ¬ 
ences, but it is always true that the group A stars are more 
centrally concentrated than the group B stars, which in turn 
are more concentrated than the group C stars. 

Table [4] shows the statistics of the radial distributions when 
the RGB stars are divided according to an empirical ridge 
line in the F343N-F555W vs. F555W diagram, instead of 
using the theoretical models as a reference. The ridge line 
(shown as a thick black line in Fig. [T]) was defined similarly 
to the ridge line for main sequence stars used for the redden¬ 
ing map, and the RGB stars were divided into two groups 
with bluer and redder colors than the ridge line, respectively. 
As in Table [3] the difference between the radial distributions 
of the two sub-samples is highly significant. For the “stan¬ 
dard analysis”, the median radii for the blue and red popula¬ 
tions are R me ;( j- = 25'.'0 + 171 and P me d,+ = 3076 ± 079 with 
P = 3.6xl0~ 4 (here, we use ’ and *+’ to denote the stars that 
are bluer and redder than the ridge line). For the brighter bins. 


Figure 9. Ratio of color distributions in the inner (R < 20") and outer (R > 
20") regions of M15. The inner regions preferentially contain stars with blue 
F343N-F555W colors. 



:l 0.104 
- 0.096 

0.088 


0.080 

0.072 

- 0.064 
| 0.056 

Figure 10. Map of the A(F343N - F555W) color of RGB stars (shown with 
black points). The mean colors are bluer near the center of Ml 5, in agreement 
with the inferred higher fraction of N-normal stars there. The dashed circle 
is centered on Ml5 and has a radius of 1'. 

the difference is again less significant because of the smaller 
number of stars, but numerically it remains consistent with the 
bins that include fainter stars and thus have better statistics. 

An alternative way of illustrating the differences in the ra¬ 
dial distributions is shown in Fig. [9] where we plot the number 
ratio of stars in the inner regions (R < 20") vs. the outer re¬ 
gions (R > 20") as a function of A(F343N - F555W). We 
see that stars with blue A(F343N - F555W) colors are more 
prevalent near the center, which is consistent with the differ¬ 
ences in the cumulative radial distributions. The difference 
between the color distributions in the inner and outer regions 
is again highly significa nt, w ith P = 6.6 x 10~ 4 . 

As discussed in Sect. |3.2| there may be small variations in 
the reddening and photometric zero-points across the field. 
Could such variations cause the observed trends? We have 
already noted that a correction for differential reddening actu- 
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Table 3 

Statistics of radial distributions for different selection criteria 


Selection 

N 

^med,A 

^med,B 

^med,C 

Fab 

Pac 

Fbc 

Standard 

+ 1 < F555W < +3 

1339 

2373 ± 172 

2974 ± l'.'O 

3274 ± 270 

1.4 x ur 3 

1.4 x 10~ 5 

0.19 

+ 1 < F555W < +2.5 

761 

2275 ± 174 

277 8 ± 175 

3171 ± 272 

0.032 

7.5 x 10~ 4 

0.23 

+ 1 < F555W < +2 

380 

2370 ± 273 

2871 ± 270 

3075 ± 2'.'3 

0.40 

0.11 

0.68 

With differential reddening correction 
+ 1 < F555W < +3 1330 22'/7 ± l'/O 

19! \ ± l'.'O 

3678 ± 274 

2.4 x ur 5 

5.1 x 10~ 9 

1.0 x 10~ 3 

+ 1 < F555W < +2.5 

760 

2270 ± 176 

2673 ± 174 

3378 ± 272 

3.8 x ur 3 

3.1 x 10~ 6 

2.2 x 10~ 3 

+ 1 < F555W < +2 

374 

2371 ± 274 

2675 ± 179 

3273 ± 272 

0.071 

5.2 x 10~ 3 

0.063 

Removing blends 
+ 1 < F555W < +3 

1272 

2571 ± 174 

3076 ± l'.'O 

3278 ± 270 

6.4 x 10“ 3 

1.1 x 10~ 4 

0.27 

+ 1 < F555W < +2.5 

730 

2372 ± 176 

2873 ± 175 

3175 ± 274 

0.081 

1.9 x 10~ 3 

0.21 

+ 1 < F555W < +2 

367 

2478 ± 276 

2875 ± 179 

3076 ± 2'.'3 

0.49 

0.19 

0.71 


Note. — “Standard analysis” refers to the default procedure described in the main text, i.e. no 
differential reddening correction or explicit removal of blends. 


Table 4 

Statistics of radial distributions split at ridge line 


Selection 

N 

^med,- 

^med,+ 

P 

Standard 

+ 1 < F555W < +3 

1364 

2570 ± 171 

3076± 079 

3.6 x 10~ 4 

+ 1 < F555W < +2.5 

769 

2377 ± 172 

2971 ± 172 

0.017 

+ 1 < F555W < +2 

381 

2573 ± 2'.'3 

2876 ± 174 

0.50 

With differential reddening correction 
+ 1 < F555W < +3 1368 23'/4 ± 079 

3079± 078 

2.7 x 10- 7 

+ 1 < F555W < +2.5 

784 

2276 ± l'.'l 

3070 ± 172 

1.5 x ur 4 

+ 1 < F555W < +2 

386 

2378 ± 177 

3071 ± 172 

0.021 

Removing blends 
+ 1 < F555W < +3 

1292 

2674 ± 173 

3176± l'.'l 

9.1 x 10~ 4 

+ 1 < F555W < +2.5 

738 

2479 ± 173 

2978 ± 172 

0.024 

+ 1 < F555W < +2 

368 

2670 ± 276 

28'.'8 ± 173 

0.46 


Note. — The columns R me:di _ and R med + give the median pro¬ 
jected distances from the center of M15 for stars to the left and 
right of the empirical ridge line, respectively. 


ally increases the statistical significance of the trends. How¬ 
ever, it is uncertain to what extent the structure in the red¬ 
dening map is real. Figure 10 shows a map of the average 
A(F343N-F555W) color across the WFC3 field. The map was 
produced in the same way as the reddening map in Fig. HI 
except that a larger smoothing radius (350 pixels) was used 
because of the smaller number of RGB stars. The mean 
A(F343N-F555W) colors are clearly bluer near the center of 
M15, which is consistent with the radial trends in Figs.[8]and 
[9] At radial distances of ~ 1' from the center, the A(F343N- 
F555W) colors are 0.02-0.03 mag redder than at the center. 
The reddest colors are seen in the lower left-hand corner, 
where Fig. [4] also shows the highest reddening, which may 
suggest that there is a real gradient in the reddening across the 
field. However, the general morphology of Fig. [To] is quite 
different from that of Fig. [4] and the former appears much 
more symmetric with respect to the center of Ml5, especially 
if one imagines subtracting an overall left-right gradient. If 
the color variations were dominated by differential reddening, 
we would expect a greater degree of similarity between the 
two maps, so it appears likely that the F343N-F555W color 
variations are, for the most part, intrinsic to M15. However, 
with our present data it remains difficult to exclude with cer¬ 
tainty the possibility that some combination of small redden¬ 
ing variations and zero-point variations could affect the ob¬ 
served radial distributions noticeably. 



Figure 11. Comparison of the SDSS A (u - g) and HST A(F343N-F555W) 
color offsets. 


3.5. Combining the HST data with SDSS photometry 

Because of the limited field of view, the WFC3 data only 
allow us to constrain radial trends out to about 130", with 
poor statistics in the outer parts because of incomplete spatial 
coverage. However, An et al. ( 2008) ) have carried out PSF- 
fitting photometry on imagesTrom the Sloan Digital Sky Sur¬ 
vey (SDSS; [York et al7[ [2000| for a number of fields around 
Galactic globular clusters, including M15. While these data 
lack the spatial resolution of HST imaging and cannot resolve 
the inner parts of the clusters, they extend to much larger radii. 
The Sloan u filter covers NH and CN absorption bands, which 
makes the u - g col or sensitive to light element abun dances. 
Th e | An et al.| ( |2008] > photometry was analyzed by Lardo et al. 
( |2011| ( who found that stars with redder it - g colors (i.e., 
the enriched stars) tended to be more radially concentrated 
in M15 (and other GCs). This trend is opposite to that seen 
in Fig. [8]from the WFC3 data in the central regions of M15. 
These results are, however, not necessarily in contradiction to 
each other, because of the limited overlap of the radial ranges 
covered by the two datasets. 

To explore the variation in the ratio of primordial vs. en¬ 
riched stars over the full radial range, we combined our HST 















































10 


S. S. Larsen et al. 


data with the SDSS photometry of|An et al.| (j2008). The 
SDSS data extend to about 2 0' from the center of Ml5, but 
following Lardo et al. ( |2011 ) we restrict our analysis to the 
innermost 10' to limit contamination by field stars. The filter 
combinations of the two datasets are similar, but not identi¬ 
cal, so we first tested whether they yield consistent results in 
the region of overlap (at radii of 60"-130"). We defined a 
A (u - g) parameter for the SDSS data in the same way as 
the A(F343N-F555W) parameter. As a reference, we used 
the same isochrone employed in Fig. [I] with the SDSS col¬ 
ors provided through the Dartmouth web i ntcr laccj^] Our 
A(m - g) parameter is closely analogous to the A„ _ g parameter 
defined by Lardo et al., the main difference being that they 
used an empirical ridge line as a reference whereas we use a 
theoretical isochrone. In Fig. 11 we compare the A(F343N- 
F555W) and A(u - g) measurements for stars in common be¬ 
tween the HST and SDSS datasets. For the SDSS data we 


adopted the 15 < g < 17 magnitude range of Lardo et al. 
(20111, where the faint limit is equivalent to F553W w 16.6 


or Mf 555 w ~ 1.24. Since this is close to the bright magni¬ 
tude cut-off in our HST data, we relaxed the bright magni¬ 
tude limit in the HST data to Mpsssw = +0.5 (for the pur¬ 
pose of this comparison only) to increase the overlap be¬ 
tween the tw o samp les, and we also eliminated the x\ cut. 
As noted by |Lardo et al.| (j201 If, the photometric errors ac¬ 
count for a significant fraction of the dispersion in u - g 
and there is indeed a substantial scatter in Fig. ITT] Never¬ 
theless, there is a significant correlation between A(F343N- 
F555W) and A(m - g), with the straight line showing a lin¬ 
ear least-squares fit to stars with -0.1 < A (u - g) < 0.1 
and -0.05 < A(F343N - F555W) < 0.25. The fit yields 
A (u-g) = (0.33 ±0.07)xA(F343N-F555W)-(0.042+0.009). 

For comparison with the fit in Fig. 11 we used 


ATLAS 12/SYNTHE synthetic spectra to calculate the expected 
color difference between N-normal and N-enhanced models 
for the SDSS filters, following the same approach as in Pa¬ 
per I and for the models in Fig. [2] For stars on the lower 
RGB, we found a difference of 0.05 mag in A (u - g), i.e. 0.31 
times the 0.16 mag difference in A(F343N-F555W). There is 
thus a very good agreement between the theoretical (0.31) and 
measured (0.33 ± 0.07) slopes of the A (u - g) vs. A(F343N- 
F555W) relation. We therefore used the linear fit in Fig. [IT 
to transform the A(m - g) values to A(F343N-F555W) values. 
To make a clear distinction between A(F343N-F555W) val¬ 
ues transformed from the SDSS photometry and those mea¬ 
sured directly in the HST data, we will denote the former by 
A(F343N-F555 W)t- in the following. 

Figure [T2]shows the number ratio Na/Nb+c as a function of 
radius forthe combined HST/WFC3 and SDSS data, where 
Na is the number of group A stars. We count the group B 
and C stars together, Nb+c , since these cannot be well dis¬ 
tinguished in the SDSS photometry, and the difference be¬ 
tween the spatial distributions of these two groups is the least 
significant according to the HST data. To account for the 
larger scatter in the SDSS data, we have extended the red and 
blue limits of the A(F343N-F555W)7- range by 0.05 mag at 
both the red and blue edge, thus including stars in the range 
-0.08 < A(F343N - F555W) r < 0.28. We have applied the 
same color cut at A(F343N - F555W)-/ = 0.05 between group 
A and B+C as in the preceding figures. Comparing with the 
photometry at larger radii, we found that contamination in the 
SDSS sample is limited to < 10% after applying these selec- 
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Figure 12. Ratio of primordial vs. enriched stars as a function of radial dis¬ 
tance. HST and SDSS datasets have been combined to cover the radial range 
from 4" to 10'. The arrow indicates the half-light radius, R/, = 1(06 iHarrisl 

Illa¬ 

tion criteria, even in the outermost bin. 

From Fig. [12] the WFC3 and SDSS samples agree reason¬ 
ably well on the N a /Nb+c ra ti° in the overlap region. At radii 
< 60", the Na/Nb+c ratio increases towards the center, con¬ 
sistent with the more centrally concentrated distribution of the 
group A stars seen in Figs. [7] and [8] Around one arcminute (3 
pc) from the center, the N a /Nb+c ratio h as a minimum and 
then increases again towards larger radii, in agreement with 
the analysis by |Lardo et al.| ([2011). A hint of this increase 
is seen already in the WFC3 photometry (although the outer¬ 
most bin has a large error bar) and from the “donot-shaped” 
morphology of Fig. [TO] but it becomes very clear in the com¬ 
bined dataset. We thus find that the Na/Nb+c ratio has a U- 
shaped dependency on radius with a minimum near the half- 
light radius (L06; Harris[l996). This is our key result. 

We note that the brightest stars in the SPSS p hotometry 
are above the threshold where Gratton et al. ( 2000j ) find evi¬ 
dence of deep mixing. These authors find modified abundance 
patterns for log L/L Q > 1.8 in metal-poor field giants, corre¬ 
sponding to M g < 0.8, whereas the absolute magnitude range 
of the M15 SDSS photometry is -0.42 < M g < 1.58. Thus 
from Figure 10 of|Gratton et al.|(J2000J), deep mixing may have 
enhanced the N abundances by up to ~ 0.5 dex in the bright¬ 
est giants in the SDSS data. This is still a relatively small 
effect compared to the overall ~ 2 dex variations between 
the different populations and, moreover, the effect would be 
to move stars from the group A to the group B+C bins and 
thereby lead us to underestimate the true Na/Nb+c ratio from 
the SDSS photometry. 


3.6. The effect of mass segregation on populations with 
different He abundance: N-body simulations 

A tendency for the enriched stars to prefer a location near 
the half-mass radius is unexpected in all current scenarios for 
the origin of multiple stellar populations in globular clusters. 
A possible explanation could be a significantly enhanced He 
abundance of the enriched stars and modification of the ini¬ 
tial radial distributions by two-body relaxation. Because of 
the He enhancement, the enriched stars would undergo faster 
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Figure 13. Observed radial distributions of N-normal stars and N-enhanced 
stars, compared with /V-body simulations for stars of different masses that 
start with the same radial distribution. The curve labeled ‘Giant stars’ is for 
masses of 0.82 M©. The lines show the simulated radial distributions after 
11.5 Gyr, normalized to the outermost measured point. 

stellar evolution so that the enriched (group B+C) post-main 
sequence stars have lower masses than the primordial (group 
A) ones for a given age. Mass segregation due to two-body 
relaxation would then push the enriched giants outwards, de¬ 
creasing their number fraction in the centre. In the outer clus¬ 
ter parts the relaxation time is much longer, preserving the 
initial ratio for a much longer time, hence one might expect to 
find a U-shaped profile in the number ratio of the populations 


similar to that seen in Fig. 12 


In order to test what mass differences between primordial 
and enriched stars are necessary to explain the radial trends 
see n in Ml5, we used the g rid o f A-body simulations used 
by [McNamara et ~ak ( |2012 i and ILiitzgend orf et al.| (2013 > 
and selected the best-fitting non-IMBH cluster from this grid. 
Fig. [13] shows the radial distribution of stars with different 
masses at T=11.5 Gyr after all clusters were scaled to the 
mass and size of M15. The “giant stars” have masses of 
0.82 M a and the lower mass stars have been arranged into bins 
of 0.50 M 0 - 0.60 M 0 and 0.60 M 0 _- 0.70 M e . Because the 
/V-body simulations of McNamara et al.|(|2012| start without 
mass segregation, we cannot expect to fit the outer profile of 
N-enhanced stars in Ml5, therefore we restrict our fit to radii 
inside the half-mass radius. Because of the long relaxation 
time in the outer regions, mass segregation due to two-body 
relaxation will not have developed there within a Hubble time, 
and a dedicated set of A'-body simulations including initial 
segregation of the different populations would be required to 
reproduce the full radial profile. 

As can be seen, the radial distribution of the giant stars in 
the A'-body simulation provides a very good fit to the distribu¬ 
tion of the N-normal giants in M15 (here, the observed radial 
distributions have been corrected for spatial incompleteness). 
The radial distribution of the most N-enhanced (group C) stars 
is best fitted by the distribution of stars with masses between 
0.50 M q < M < 0.60 Mq in the A'-body simulation, whereas 
the intermediate stars (group B) correspond better to masses 


of 0.60 M q < M < 0.70 Mq. Hence if the different radial dis¬ 
tributions are due to two-body relaxation driven mass segre¬ 
gation, then the most enriched RGB stars would need to have 
masses of about 0.25 M 0 less than the primordial giants. This 
requires an extreme degree of helium enrichment: for exam¬ 
ple, according to the models of Dotter et al.](|2007(l, Y = 0.40 
instead of Y - 0.25 produces a mass difference of 0.20 M G 
for stars on the lower RGB at an age of 13 Gyr (M = 0.79 M 0 
for Y = 0.25 and M - 0.59 M 0 for Y = 0.40). Because the 
initial conditions of the simulations are not mass segregated, 
the required mass differences (and, consequently, the required 
He enhancement) may be considered conservative estimates. 
If the cluster started out with a more concentrated enriched 
population, an even larger difference in mass between the gi¬ 
ants of the two populations would be required to reverse this 
and produce the observed differences in the available time. 

3.7. Observational constraints on He abundance variations 

in Ml 5 

A significant enhancement of the He abundance in the 
enriched stars is expected in most scenarios, due to the 
H-burning nucleosynthesis involved ([Ventura et al.l|2001| 


D’Antona et al.l2002HDecressin et al.|2007b|lde Mink et al.| 

2009|l. The observations of distinct, parallel main sequences 
that cannot be explained by differences in heavy element 
abundances provide strong evidence for He abundance vari¬ 
ations in some GCs. In some clusters, such as a> Cen and 
NGC 2808, the He fraction of the enriched stars may be as 
high a s Y ~ 0.40 ( Bedin et al.|2004 Norris|2004 Piotto et al. 
2007), but in others the enhancement is much more modest; 
photometry of the clusters 47 Tuc, NGC 6397, and NGC 6752 
indicates He enhancement in the range AT = 0.01-0.03 
( |Milone et al.||2012a|b[ |2013] >. It is difficult to measure He 
abundances directly from spectroscopy of cool stars because 
of the large difference between the two lowest energy levels 
of neutral He. In hot horizontal branch stars, where He is 
more readily measurable, the surface composition may have 
been heavily modified by stellar evolutionary effects (Behr 
|et al.j[2000{ Valcarce et al.||2014| . Star-to-star differences in 
the strength of the chromospheric He i 10 830 A line in red 
giants have been observed in a> Cen and NGC 2808 and im¬ 
ply He abundances consisten t with the large variations de¬ 
rived from photometry (Dupree et al. |2011| |Pasquini et al. 


201 1 ), but deriving accurate abundances from the line is difii- 
cult. Unfortunately, neither high-precision photometry of suf¬ 
ficient depth to reveal multiple main sequences nor spectro¬ 
scopic constraints on the He abundance of red giants currently 
exist for M15. 


As noted in Sec. 3.1 not only the main sequence, but 
also the effective temperatures of stars on the lower part of 
the RGB should be sensitive to He abundance (Salaris et al. 


[2006] [Pi Criscienzo et al. !20l 1 Beccari et al. 2013| . Fig. 14 

shows [Potter et al.| ( |2007[) isochrones for an age of 13 Gyr 
and [ayFe] = 0.4 and different He abundances (T = 0.25, 
0.33, and 0.40). Here we have scaled [Fe/H] for T = 0.33 
and T = 0.40 to the same total Z (i.e., [Fe/H] = -2.25 
for T = 0.33 and [Fe/H] = -2.2 for T = 0.40). In both 
F555W-F814W and F343N-F555W, the colors of RGB stars 
are predicted to be bluer for the He-enhanced models (keep¬ 
ing the light element mixture fixed). At M F 555 -w = +2, the 
difference between the T = 0.25 and T = 0.40 models is 
-0.036 mag in F555W-F814W and -0.034 mag in F343N- 
F555W (for T = 0.33 vs. T = 0.25, the corresponding dif¬ 
ferences are -0.020 mag in F555W-F814W and -0.021 mag 
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Figure 14. |Dotter et al. 1(2007] ! isochrones for different He abundances. The 
isochrones have an age of 13 Gyr and the same Z, corresponding to [Fe/H] = 
-2.3, [Fe/H] = -2.25, and [Fe/H] = -2.20 for Y = 0.25, Y = 0.33, and 
Y = 0.40, respectively. 


in F343N-F555W). If we instead keep the iron abundance rel¬ 
ative to hydrogen fixed at [Fe/H] = -2.3, then the F555W- 
F814W color offsets remain the same as above, whereas the 
difference between the Y = 0.25 and Y = 0.40 colors now 
amounts to -0.050 mag in F343N - F555W. In any case, the 
observed spread in F343N-F555W cannot be explained by He 
abundance variations and remains dominated by N abundance 
variations (cf. Fig. [2}. For F555W-F814W, however, the sit¬ 
uation is the opposite: this color is largely insensitive to N 
abundance and any difference in F555W-F814W would there¬ 
fore be attributable to He abundance variations. 

In Fig. fl5|we show the observed A(F555W-F814W) dis¬ 
tributions for the group A, B, and C stars. The mean col¬ 
ors of the three groups are <A(F555W - F814W)) A = 0.015 
mag, (A(F555W - F814W)) B = 0.011 mag, and <A(F555W - 
F814W))c = 0.007 mag. There is, indeed, a tendency for 
the group B and C stars to have bluer F555W-F814W col¬ 
ors than the group A stars, as would be expected if there is a 
difference in He abundance. From Fig.fl4l we get A(F555W- 
F814W)/AT = 0.24 (at = 2), so a color difference of 

0.008 mag between the group A and C stars corresponds to 




Figure 16. Distributions of A(F606W-F814W) for group A, B, and C stars. 
There are no significant differences between the three groups. 


AT ~ 0.033. The corresponding mass difference on the RGB 
is only 0.045 M Q , which is much too small to explain the dif¬ 
ferences in spatial distribution as a consequence of mass seg¬ 
regation. The differences in (A(F555W - F814W)) between 
the groups become even smaller (< 0.002 mag) if the differ¬ 
ential reddening correction is applied. 

As an independent check, we made the same compari¬ 
son using the ACSGCS F606W-F814W colors. Unlike the 
F555W-F814 colors, this filter combination is fully indepen¬ 
dent of the F343N-F555W colors. It does have a slightly nar¬ 
rower color baseline than F555W-F814W and the predicted 
color difference between T = 0.25 and T = 0.40 RGB stars 
is reduced to -0.028 mag at M P g 06 w = +2 (the T = 0.25 vs. 
T = 0.33 difference is -0.015 mag). Nevertheless, this should 
still be easily detectable. We matched our WFC3 data with the 
ACSGCS photometry and defined a A(F606W-F814W) index 


in the same way as for the WFC3 observations. Figure 16 
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shows the distributions of A(F606W-F814W) colors for the 
group A, B, and C stars. We find <A(F606W - F814W)) A = 
-0.013 mag, <A(F606W - F814W)) B = -0.009 mag, and 
(A(F606W - F814W))c = -0.012 mag. Note that the group 
A and C stars now have very similar F606W - F814W colors, 
and the group B stars differ by only 0.004 mag from the group 
A stars. Formally, this corresponds to a difference in the He 
abundance of AT ~ 0.021, but the statistical significance of 
any differences between the color distributions of the three 
groups is small: a K-S test yields Pab = 0.05, Pac = 0.60, 
and P BC = 0.15, respectively. 

From these comparisons, we conclude that a conservative 
upper limit on the He abundance variations in M15 is AY < 
0.03, which corresponds to a mass difference of < 0.04 M 0 
on the RGB for coeval populations. As recently emphasized 
by Dotter et al. ( j2015l> , absolute constraints on He abun¬ 
dance variations from CMD analyses remain uncertain. How¬ 
ever, even with a generous allowance for model uncertain¬ 
ties, it seems difficult to accommodate the large variations in 
He abundance that are necessary in order for dynamical ef¬ 
fects to produce the observed radial distributions of the sub¬ 
populations in M15. 

4. DISCUSSION 

The radial trends in the ratios of N-normal to N-enhanced 
stars found in the preceding sections are very surprising, given 
that all current scenarios for the origin of multiple stellar pop¬ 
ulations predict that the enriched population should be more 
centrally concentrated. While the exact radial distributions 
predicted by each of the models are clearly subject to consid¬ 
erable uncertainty, it is difficult to imagine formation scenar¬ 
ios in which the enriched stars preferentially avoid the center, 
as appears to be the case in M15. 

In the absence of an obvious explanation related to the for¬ 
mation of the populations, we have considered the possibility 
that dynamical evolution is responsible. We have argued that 
mass segregation might produce trends similar to those ob¬ 
served in the inner regions, but only if the enriched stars are 
very strongly He-enhanced. At 13 Gyr, the difference between 
models for giants with normal (Y = 0.25) and the most He- 
enhanced composition available (Y = 0.40) is about 0.2 M 0 , 
which is barely sufficient - our A'-bodv simulations indicate 
that a mass difference of ~ 0.25 M Q or more is required. Fur¬ 
ther out, any trends set up at formation would be preserved. 
Unfortunately, the currently available constraints on the He 
abundances of the different populations do not appear to sup¬ 
port this explanation. Instead, we estimate that the masses 
of stars on the lower RGB differ by less than 0.04 M 0 , effec¬ 
tively ruling out mass segregation via two-body relaxation as 
a viable explanation for the observed radial trends. 

Differences in the masses of red giants could also be pro¬ 
duced by age differences, but in order for the enriched giants 
to be pushed outwards by mass segregation, they would have 
to be older than the primordial ones. Furthermore, the age dif¬ 
ference would have to be extremely large; even a difference of 
3 Gyr corresponds to a mass difference of only 0.05 M 0 at the 
turn-off. Exploring this possibility fully would probably re¬ 
quire a dedicated A'-body simulation that explicitly takes the 
age differences into account, but large age differences appear 
to be ruled out by the narrow subgiant branch in M15. 

4.1. Other constraints on He abundance 

It should be kept in mind that the relations between He 
abundance and the broad-band colors on the RGB remain to 


be verified observationally, although they appear to be a ref 


atively solid prediction of the models. Di Criscienzo et al. 
( |2011[ ) found color differences between He-normal an d He^ 
enhanced RGB stars similar to those predicted by th e |Dot-| 
|ter et al-1 ( |2007| ) mod els, as did |Beccari et ah] ( |2013| l based 
on BaSTI isochrones (Pietrin ferni et al.|2004) . From a self- 
consistent modeling of isochrones with enhanced He (T = 


+0.35) and modified light element abundances, Salaris et al. 
(2006[) found a difference in V - / of ~ 0.03 mag with respect 


to He-normal stars on the lower part of the RGB, again similar 


to the offsets in Fig. 14 


It would certainly seem worthwhile to obtain deep, high 
precision photometry of the main sequence in M15 in or¬ 
der to obtain better constraints on any variation in He abun¬ 
dances. In the meantime, another indicator worth exploring is 
the horizontal branch (HB) morphology. The suggestion that 
HB morphology and He abundance of GCs may be related 
is as old as the second parameter problem itself (Sandage & 


Wildey 

1967||van den Bergh 1967| Kraft||1979| Freeman & 

Norrisf] 

with an 

981). M15 has a complex, bimodal HB morphology 

extended blue tail that has proven difficult to model in 

detail (Buonanno et al. 19851 Du ire II & Harris 1993 [Moehler 


et al.|1995 1 . From the HB morphology and the RR Lyrae pe 

riod distribution, [D’Antona & C aloi (2008) inferred a second 
generation fraction of 80% and a moder ate He enhancemen t 
(Y = 0.26-0.30). Using the isochrones of Dotter et al. |(j2007J), 
an Y — 0.30 RGB star with an age of 13 Gyr has a mass of 
about 0.72 M 0 , which is a difference of -0.07 M 0 with r espect 
to the standard (Y = 0.25) composition. Jang et ak|(|2014i 
suggested that the stars in the blue tail of the M15 HB belong 
to a population with Y = 0.33 that accounts for 42% of the 
stars in the cluster, but their models assume a 1 Gyr age dif¬ 
ference between the first and subsequent generations of stars. 
In this case, the mass of a He-rich RGB star would be 0.69 
M 0 . From Figure [13] even a ~ 0.1 M 0 mass difference be¬ 
tween the different populations seems insufficient to explain 
the observed radial trends as an effect of mass segregation. 

While recent work has focused on He, it is well known that 
there are other parameters, such as the abundances of CNO, 
that can affect HB morphology flHartwick & McClure||1972| 


Dorman et al.|199ll|Dorman|19 92; Salar is et al.|20 06; Milone 

et al.||2014| i. Indeed, |D’Antona & Caloi| ( |2008| l found that 
they were not able to get a satisfactory fit to the M15 HB 
by varying only He. Data that allow comparison of He abun¬ 
dances derived from the HB and other methods are only avail¬ 
able in a few cases. In 47 Tuc, |D’Antona & Caloi] ( |2008[ > 
found a second-generation fraction of ~ 25% a nd Y — 0.27- 
0.32 (AT = 0.02-0.07) from the HB, whereas [Milone et al.| 
(j2012b) found AT ~ 0.015 and a second-generation fraction 
of ~ 70% from analysis of the full CMD, thus favoring a He 
enhancement toward the lower end of the range indicated by 
the HB analysis. Similarly, Di Criscienzo et al.j( j2010| > found 
that the HB morphology of 47 Tuc and spread in the luminos¬ 
ity of stars on the sub-giant branch could be explained by a 
small but real He abundance spread of AT = 0.02, in agree¬ 
ment with the analysis of Milone et al. (j20T2b|, combined 
with an enhancement of the C+N+O sum in a fraction of the 
He-enhanced stars. In NGC 6397, [DAAntona & Caloi| (2008) 
suggested a negligible first generation fraction and T = 0,28 


(AT = 0.04) for the second generation, whereas Milone et al. 
(j2012ajl found that 70% of the stars belong to an enriched 
population with AT ~ 0.01 (T ~ 0.26). Finally, in the case 
of NGC 2808, |Dalessandro et al.j ( |201 1[ > found that the HB 
morphology could mostly be well matched by a model that 
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Figure 17. Ratio of primordial to enriched stars as a function of proje cted 
dist ance for the best-fitting A-body model to M15 as described in Sect. |3.6| 
and |4.2| For this plot we assumed that the stellar evolution of both stellar 
populations is the same. All stars with initially low and high orbital ener¬ 
gies were assumed to be primordial stars while stars with intermediate orbital 
energies were assumed to be enriched stars. Despite orbital mixing due to 
two-body relaxation, a significant difference between the two populations is 
preserved until the end of the simulation at T=11.5 Gyr. 

incorporates the constraints on He abundance from the three 
distinct main sequences in this cluster, although this model 
had some difficulty reproducing the hottest “blue hook” stars 
on the HB. These comparisons underline the considerable un¬ 
certainties involved in inferring He abundances from CMD 
analyses, and from the HB in particular, but also show that 
the extreme degrees of He enhancement in clusters like a> Cen 
and NGC 2808 are far from universal. 


4.2. Mixing of sub-populations 

Any initial segregation of equal-mass sub-populations 
within a GC is expected to be eventually erased by mixing 
due to two-body relaxation (|Decressin et al.|2008). However, 
this is a relatively slo w process; the A'-body simulations of 
Vesperini et al. ( |2013| ) indicate that any initial differences in 
the half-mass radii of different populations should remain de¬ 
tectable at least until the cluster is 10 half-mass relaxation 
times old (here referring to the current half-mass relaxation 
time). The current half-mass relaxation time of M15 is about 
10 9 years (Djorgovski 1993) so it does not seem unreasonable 
that the relatively moderate variations in the number ratios (by 
a factor of 2-3) observed in the central regions of M15 are still 
preserved. 

In order to check how quickly initial segregation is erased 
by two-body relaxation if the two sub-populations have equal 
mass, we agai n used the best-fitting model from the grid of 
|McNamara et al. ( j2012j ) and assumed that all stars with ini¬ 
tial orbital energies between -1.0 < E < -0.2 (in A'-body 
units) are enriched stars while all stars with lower and higher 
energies are primordial stars. Figure ITT] shows the projected 
number ratio of M > 0.6 M Q primordialto enriched stars that 
corresponds to this choice at the start of the simulations and 
after 12 Gyr of evolution. Due to two-body relaxation and re¬ 


sulting orbital mixing, the initial segregation between the two 
components is strongly weakened but not completely erased 
and the number ratio between the populations agrees qualita¬ 
tively with what we find for M15 as shown in Fig. [12] We find 
similar amounts of remaining segregation for otherchoices of 
initial segregation, which demonstrates that the signal that we 
see is not due to one particular choice of initial segregation. 
The presence of a segregation of different stellar generations 
in M15 is therefore not in conflict with the idea that the pop¬ 
ulations have the same mass. 


4.3. Are there other cases like Ml5? 

While the enriched population has generally been found to 
be more centrally concentrated than the primordial population 
within GCs, the majority of the studies to date have not looked 
at the central regions of clusters, largely due to resolution ef- 


fectsr For example, Carretta et al. 
Kravtsov et al.[(2011[), Johnson & 


([20T0|, ILardo et al.| <f20Tl|), 
fe Pilachowski (2012), and 


Beccari et al.|(|2()l 3| have all used ground based observations 
(photometry and/or spectroscopy) to study the relative spatial 
distributions of the populations, and have concluded that the 
enriched population is significantly more centrally concen¬ 
trated than the primordial one. The above studies have been 
focused outside the central ~ 1' for their respective clusters, 
typically corresponding to 1 - 2 pc. In Paper I we found hints 
of the enriched population being more centrally concentrated 
than the primordial one in the four metal-poor GCs in the For¬ 
nax dSph, but again this was restricted to radii outside 1/2-1 
half-light radii. As seen in our analysis of M15 (also studied 
by jLardo et al.|2011| >, outside this radius the enriched popu¬ 
lation is indeed more centrally concentrated, however inside 
this radius a reversal occurs. Hence, it is possible that such 
a reversal in the enriched/primordial population ratios in the 
central regions is a relatively common feature, and has gone 
undetected due to resolution constraints. 

A few HST studies have found that the enriched popula¬ 
tion does remain more centrally concentrated even in the core 
of the cluster. Bellini et al.| ((2009)) used HST imaging to 
study the central regions ofw Cen, and found that the en¬ 
riched population does (slightly) dominate in the inner two 
core radii, but outside this radius, the enriched population is 
significantly more centrally concentrated that the primordial 
one. Milone et al. (j2012b| found that the enriched/primordial 
ratio increases toward the center in the massive GC 47 Tuc, 
although they have only two bins within the half-mass radius. 

Finally, we note that other cases of more centrally concen¬ 
trated primordial populations may already have been found. 
In NGC 2419, |Beccari et al.| (|2013J) found that giants with 
blue u - V colors are more centrally concentrated than those 
with redder u — V colors, i.e., the a - V colors in NGC 2419 
show the same behavior as the F343N-F555W colors in M15. 
Beccari et al. attributed the color differences to He abundance 
variations (so the blue stars would correspond to the enriched 
population), but also found that stars with anomalous (i.e., 
depleted) Mg abundances tended to have redder than average 
u-V colors. It would seem, therefore, that an alternative inter¬ 
pretation of their observations is that the giants with blue u-V 
colors are, in fact, stars with normal (primordial) composition, 
and stars with red u-V colors are enriched stars. NGC 2419 
would then be similar to M15 in having a more centrally con¬ 
centrated primordial population. However, it should be noted 
that the behavior of Mg in NGC 2419 is somewhat unusual, 
with [Mg/Fe] reaching very low values in some stars, and 
no obvious (anti-)correlation between [Na/Fe] and [Mg/Fe]. 
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ing a large variation in [K/Fe] 

Cohen & Kirby 2012, Ven- 

tura et al. 2012{ 

Beccari et al.|2( 

313). Although the enriched 

stars in NGC 24 

19 might be more He-enhanced than those 


in M15, NGC 2419 is among the most extended GCs in the 
Milky Way and has very long central and half -mass relaxation 
times (t rc ~ 10.5 Gyr and r,y, ~ 19 Gyr; |Djorgovski|[T993|l. 
It shows no evidence for significant mass-segregation (Da- 


lessandro et al. |2008[|Baumgardt et al.|2009 ;, Bcllazzini et al. 

2012 | >, and it might therefore also be difficult to explain dif- 

ferences in the radial distributions of sub-populations within 
NGC 2419 dynamically even if the populations have signifi¬ 
cantly different He abundances. 

A recent ground-based study by Alonso-Garcfa et al. (in 
prep.), using Stromgren photometry of NGC 288, also found 
that the primordial population is more centrally concentrated 
than the enriched population, from 015 — 7' (~ 1.3-18 pc) 
from the cluster core. 

These cases, along with our results for M15, show that GCs 
display a wide variety of behaviors in the relative distributions 
of the enriched and primordial stars, in apparent contradiction 
with the standard scenarios for the origin of multiple popula¬ 
tions within GCs. 


5. SUMMARY AND CONCLUSIONS 

We have combined HST/WFC3 F343,F555W,F814W and 
SDSS u,g observations to study the radial distributions of red 
giants with N-normal and N-enhanced composition over a ra¬ 
dial range of 4"-600" (~ 0.06 - 10 half-light radii) in the 
globular cluster M15. Our findings are as follows: 

• The spread in the F343N-F555W colors of RGB stars of 
a given magnitude is far greater than the observational 
errors and implies a variation in [N/Fe] of about 2 dex, 
which is consistent with previous spectroscopic results. 

• Dividing the stars into three groups with “primor¬ 
dial” (i.e., similar to halo field stars), intermediate and 
strongly enriched nitrogen abundances (group A, B, 
and C), we find that the group A stars are the most 
centrally concentrated within the WFC3 field of view 
and the group C stars the least centrally concentrated. 
The difference is highly significant and contrary to the 
expectations from current scenarios for GC formation, 
which predict that the stars with primordial composi¬ 
tion should be the least centrally concentrated. The 
group B stars have a degree of central concentration in¬ 
termediate between the A and C groups, but the differ¬ 
ence between the B and C stars is less significant. 


• However, the small differences in optical colors on the 
lower RGB suggest that there are no large differences in 
the He abundances of primordial and N-enriched stars 
(AT < 0.03), with a corresponding mass difference of 
less than 0.04 M 0 if the stars have the same ages. 

We are thus left with no convincing explanation for the 
observed radial distributions of different stellar populations 
in M15. We find no evidence that variations in the fore¬ 
ground reddening might cause the observed trends; a differ¬ 
ential reddening correction based on the F606W-F814W col¬ 
ors of main sequence stars actually increases the significance 
of the trends. However, small color variations across the field 
might also be caused by instmmental effects. Data in more 
passbands would be required in order to quantify these effects 
better. 

If the overall trends found here, including the lack of a sig¬ 
nificant difference in He abundance, are confirmed (for ex¬ 
ample by deep, multi-passband photometry of the main se¬ 
quence), then the differences in the central regions are un¬ 
likely to be of a dynamical origin and, presumably, must re¬ 
flect the conditions at the time of formation. We have shown 
that such differences could be preserved until the present 
epoch, due to the relatively slow nature of orbital mixing. 
This would represent a challenge to all current scenarios for 
the formation of GCs. 
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